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\ Recent investigations on the central stars of planetary nebulae (CSPN) indi- 

cate that the masses based on model atmospheres can be much larger than the 
. masses derived from theoretical mass-luminosity relations. Also, the disper- 

tin ■ sion in the relation between the modified wind momentum and the luminosity 

' depends on the mass spread of the CSPN, and is larger than observed in 

massive hot stars. Since the wind characteristics probably depend on the 
metallicity, we analyze the effects on the modified wind momentum by con- 
sidering the dispersion in this quantity caused by the stellar metallicity. Our 
CSPN masses are based on a relation between the core mass and the nebular 



abundances. We conclude that these masses agree with the known mass distri- 
■ bution both for CSPN and white dwarfs, and that the spread in the modified 

wind momentum can be explained by the observed metallicity variations. 

Key Words: Planetary nebulae — Stars: Planetary nebula central stars — 
Stars: Mass loss 
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CD . 1. INTRODUCTION 

co 

, Stellar winds are observed in practically all the HR diagram, from dwarf, 

solar-like stars to the hot supergiants and massive stars in the main sequence. 
In particular, AGB stars have slow, massive winds leading to the ejection of 
the stellar outer layers as a planetary nebula (PN). At this stage, the stellar 
remnant is extremely hot, originating a fast, massive wind which is similar to 
. the radiative winds observed in hot, massive young stars. These winds are 

observed in central stars of planetary nebulae (CSPN), both for H-rich stars 
and for stars with Wolf-Rayet characteristics. The origin of the winds of these 
stars is attributed to the stellar radiation pressure on metallic lines, a similar 
mechanism known to be active in hot, massive galactic stars. 

In the last few years, some work has been done in order to use the available 
information on the stellar winds of CSPN to determine physical properties 
of the stars, such as th eir radii, luminos ities a nd masses (see for example 
Kudritzki & Puis bood Pauldrach et al. l2003h . In particular, observed or 



inferred properties such as the mass loss rate and the wind terminal velocity 
can be used in order to constrain some physical properties of the stars, such as 
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their luminosities and masses. As a result, there seems to exist a discrepancy 
between recent mass determinations for CSPN based on model atmospheres 
affected by winds and the corresponding determinations from mass-luminosity 
relations derived from evolutionary tracks of post-AGB stars. There is a much 
larger mass spread obtained from model atmospheres as compared with the 
predictions of standard mass-luminosity relations and with the known mass 
distribution of CSPN and white dwarfs. 

For a stellar wind with mass loss rate M — dM/dt and terminal velocity 
v oo in a star with radius i?» , a correlation can be obtained between the mod- 
ified wind momentum, defined as p w = M Voo y/R* and the stellar luminosity 
L/Lq. This correlation is well defined for massive stars, and is also approxi- 
mately valid for CSPN. However, use of a standard mass-luminosity relation 
implies a relatively large spread in the modified wind momentum at a given 
luminosity for CSPN, in contrast with the results from model atmospheres, 
which show a be tter a greement with th e rela tion derived from massive stars 



(Kudritzki et al. 119971 ; Pauldrach et al. 120041 ). 

In the present work, we analyze the available information on the wind 
momentum in a sample of well-studied CSPN in the Galaxy, and use this 
information in order to constrain basic stellar properties. In particular, we 
look for correlations involving the observed wind and nebular properties, es- 
pecially the chemical abundances. In fact, the stellar metallicity (and as a 
consequence the nebular chemical composition) has some effect on the wind 
properties, which can be seen for example by comparing galactic stellar winds 
with the corresponding quantities in stars of the Magellanic Clouds, which are 
more metal-poor than the Milky Way. This can be explained by the fact that 
the radiative mechanism operates essentially on the metal absorption lines, 
so that the stellar metallicity plays a role in the mass loss process. On the 
other hand, there is a clear correlation between the nebular chemical com- 
position and the stellar mass, especially regarding those element ratios that 
are affected by the evolution of th e PN progenitor stars, such as the N/O 



and He/H ratios (cf. Maciel 120001 ) . As a consequence, the observed wind 
properties, such as the wind momentum or some related quantity should in 
principle be related to stellar mass and metallicity. The IAG/USP group has 
a considerable experience in the determination and analysis of the chemical 
composition of PN in the Galaxy and the Magellanic Clouds (see for example 



Costa et al. 120041 ) and references therein. In this work we use this database 
and examine the effects of the nebular metallicity in the observed properties 
of the corresponding stellar winds in a sample of well studied galactic CSPN. 

2. MASSES OF CSPN 
2.1. Previous mass determinations 



As recently discussed by Pauldrach et al. (|2003l . 120041 ) . there is a clear 
discrepancy between the mass-luminosity relation for CSPN as derived from 
their improved model atmospheres and the relation previously obtained by 
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Fig. 1 . Luminosities and masses of CSPN. E mpty c ircles: data by Kudritzki et al. 
(1997), iilled circles: data by Pauldrach et al. |2004h ; stars: N/O masses (see text). 



Kudritzki et al. (|l997j ). which was obtained from the application of post-AGB 
evolutionary tracks. The new results produce a larger spread in the stellar 
masses, 0.4 < M/M Q < 1.4 , as compared with the range 0.6 < M/M Q < 1 
from Kudritzki et al. (|1997l ). Moreover, the predicted luminosities are lower 
for a given mass in the new models, and the overall correlation is less well 
defined, as can be seen fro m Fig . 1. In this figure, empty circles are the 
results by Kudrit zki et al. <| 19971) and the filled circles refer to the data by 
Pauldrach et al. |2004l ). The objects included in the figure are NGC 2392, 
NGC 3242, IC 4637, IC 4593, He2-108, IC 418, Tc 1, He2-131, and NGC 6826. 
The values of the corresponding quantities (stellar radius, mass, luminosity, 
terminal velocit y, mas s loss rate, etc.) are given in the original papers by 
Kudritzki et al. |l997t ) and Pauld rach et al. (20031 : |2004) . 

The work of Pauldrach et al. (|2003l ; 120041) is based on hydrodynamically 
consistent, spherically symmetric model atmospheres which are able to repro- 
duce the observed ultraviolet spectra of hot, massive stars. On the basis of 
the observed wind characteristics, namely the terminal velocity and the mass 
loss rate, the main stellar parameters can be derived. In this case, the ob- 
tained stellar mass is independent of any adopted mass-luminosity relation, 
which in principl e may be an advantage. On the other hand, the results by 
Kudritzki et al. (1997) are based on a core mass-luminosity relationship for 
post-AGB stars, in the framework of a modelling of stellar H and He line 
profiles. Although this method may be affected by complications such as a 
contamination by nebular lines, it generally provides accurate determinations 
of the wind properties as well as of the stellar gravity and mass. Therefore, 
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we need further information in order to clarify the discrepancy in the stellar 
masses observed in Fig. 1. 



2.2. The N/O masses 

This problem can be investigated taking into account the CSPN masses 
derived fr om ob served nebul ar abu ndances of N/O, as discussed in Cazetta 
& Maciel <|2000h and Maciel (|2000l ). From AGB evolution theory, there is a 
general correlation between the N/O abundan ce rat io and the core mass, that 
is, the CSPN mass (see for example Marigo 120011 ). Theoretical models for 
AGB stars predict that higher stellar masses are associated with larger N/O 
abundances, which are enhanced by the dredge-up e pisode s that occur in the 
CSPN, especially the second one. Cazetta & Maciel (2000) and Maciel (2000) 
presented a detailed discussion of the core mass-N/O abundance relation based 
both on theoretical models of AGB stars and on empirical mass determina- 
tions. This relation takes into account several initial mass- final mass relations 
available in the literature, as well as different individual mass determinations, 
so that it can be considered as independent of a sing le m ass-luminosity rela- 
tion, such as the one adopted by Kudritzki et al. (|1997T ). As a result, two 
calibrations were determined, which were referred to as the high-mass and 
the low-mass calibration (cf. Maciel |2000| for details). The latter is consid- 
ered as more accurate, as it produces core masses M» > 0.55 M and main 
sequence masses Mmfs > 1 M©, in agreement with the detailed CSPN masses 
of Stasinska e t al. (| 19971) as well as with the masses of Type II PN proposed 
by Peimbert ( 1978h . The adopted initial mass-fin al mas s relation was based 
on the gravity distance work of Maciel & Cazetta (|1997l ). 

It is then interesting to investigate whether the masses obtained from the 
analysis of nebular abundances can be used to distinguish between the masses 
derived from model atmospheres and s tellar wind a nalysi s. 

Following Cazetta & Maciel (|2000h and Maciel (|2000l ). the core mass (in 
solar masses) can be written as 



M* = 0.7242 + 0.1742 log (N/O) 
for -1.2 < log (N/O) < -0.26 and 



(1) 



M* = 0.825 + 0.936 log (N/O) + 1.439 [log (N/O)] 2 , (2) 

for —0.26 < log (N/O) < 0.20, where (N/O) is the nitrogen over oxygen abun- 
dance ratio by number of atoms. Columns 3 to 6 of Table 1 show the abun- 
dances of oxygen, nitrogen, and N/O adopted from our IAG/USP database, 
and the calculated masses for the sample of CSPN shown in Fig. 1. The object 
IC 4637 is not in our database, so that we have adopted an average abundance 
log(N/0) = —0.56, which is typic al for elliptical planetary nebulae, according 
to the discussion of Gorny et al. (|1997l ). 



Fig. 2 shows the masse s as given by Pauldrach et al. (|2004 solid dots) 
and Kudritzki et al. (|1997l empty circles) as a function of the masses given in 
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Name 


e(O) 




loefN/01 




l0f?"f Tj* 1 Jjr^\ 


lot? T)»„ 


1 


NGC 2392 


8.50 


8.38 


-0.12 


0.73 


4.08 


26.63 


2 


NGC 3242 


8.66 


7.91 


-0.75 


0.59 


3.85 


26.44 


3 


IC 4637 






-0.56 


0.63 


3.91 




4 


IC 4593 


8.32 


7.24 


-1.08 


0.54 


3.77 


25.98 


5 


He2-108 


8.22 


7.85 


-0.37 


0.66 


3.96 


26.17 


G 


IC 418 


8.54 


7.82 


-0.72 


0.60 


3.87 


26.35 


7 


Tc 1 


8.71 


7.67 


-1.04 


0.54 


3.77 


26.37 


8 


He2-131 


8.67 


8.23 


-0.44 


0.65 


3.95 


26.60 


9 


NGC 6826 


8.43 


7.50 


-0.93 


0.56 


3.80 


26.14 



Table 1, which we may call N/0 masses. It can be seen that the N/0 m asses 
are generally clo ser to the masses derived by Kudritzki et al. 1 1997T) and 
Kudritzki et al. d2006h than to the values by Pauldrach et al. (|2004h . which 
show a much higher mass dispersion. In view of the similarity betwee n our 
N/O masses and the most of the masses obtained by Kudritzki et al. (|1997 , 
it is tempting to calibrate the luminosities using the mass-luminosity relation 
used in that work, which can be approximately written as 



log(L„/L Q ) ~ 2.90 + 1.61 (M„/M Q ) 



(3) 



The derived luminosities are given in column 7 of Table 1 and the N/O masses 
and corresponding luminosities are also included in Fig. 1 as black stars. 
It can be seen that these masses have a much more limited range, roughly 
0.5 < M(Mq) < 0.8, than obtained by Pauldrach et al. (120041) . and are also 
somewhat more restricted than the values by Kudritzki et al. (|1997l ). 



2.3. Discussion 



A comparison of the CSPN masses discussed in the previous subsections 
can be made both on the basis of the known mass distributions of these objects 
and their white dwarf descendants as well as of the rather limited individual 
mass determinations for the central stars. 



A set of central star masses was obtained by Gorny et al. (jl997l ). based 
on measured nebular expansion velocities and central star parameters, within 
the framework of a simple evolutionary model for planetary nebulae. A com- 
parison of our derived masses given in Table 1 with their results shows a good 
agreement, as all objects have a mass difference under 0.10 Mq, and five of 
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Fig. 2. Masses of CSPN asa function of the N/O masses give n in Table 1. Empty 
circles: Kudritzki et al. (|1997T ), filled circles: Pauldrach et al. l| 20041 ). 



the stars have even smaller differences, under 0.03 M©. The only exception is 
NGC 2392, for which our derived mass is 0.12M© larger than the mass given 
by them, a difference of about 17%. However, this object has a very high 
nitrogen (and He, see Section 3) ab undan ce, which is consistent with a larger 
mass than derived by Gor ny et a l. (|1997 ). 

Recently, Zijlstra et al. (2008) presented the results of a 25-year monitoring 
program of radio observations of the planetary nebula NGC 7027, possibly 
the best studied object of its kind. From the evolution of the radio flux 
during this period, an improved expansion distance was obtained. Taking 
into account theoretical tracks for post-AGB stars, a distance-independent 
mass of 0.655 ± 0.01 M© was derived. This result agrees very well with our 
N/O mass for this object, which is 0.66 M© , ado pting abundances fr om ou r 
database (see for example Maciel & Quireza 1999, Maciel & Chiappini ll994T ). 
Using alternative models, Zijlstra et al. (2008) find a slightly larger mass, still 
in good agreement with our result. 

Also, in a recent work Traulsen et al. (J2005J) analyzed high resolution ul- 
traviolet spectra of a few hydrogen-rich CSPN obtained with HST and FUSE, 
and were able to locate these objects in a logg x logT e // plane using theo- 
retical evolutionary tracks. For the three objects in our abundan ce dat aset 



2004 and 



(NGC 4361, NGC 6853, and NGC 7293, see also Perinotto et al. 
Henry et al. 2004 ) , the central star masses obtained by Traulsen et al. ( 20051 ) 
cluster around 0.6 M©, within about 15% of the N/O masses derived for these 
objects. 
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Stasinska et al. (|1997l ) applied the method described by Gorny et al. 



(|1997h and derived a mass distribution of CSPN with a very restricted range, 
in which more than 80% of the objects have masses between 0.55 and 0.65 Mq. 
This distribution peaks around 0.60 Mq with some dispersion that depends on 
the adopted nebular mass, and falls steeply towards large masses. Therefore, 
these results fit well the lower masses found by the N/ O meth od, so that the 
large masses found in a few cases by Kudritzki et al. (|1997h and especially 
by Pauldrach et al. (|2004l ) seem rat her u nlikely. A similar conclusion has 
been recently reached by Napiwotzki (|2006r ). who argued that the high CSPN 
masses close to the Chandrasekhar limit derived by Pauldrach et al. (|2004r ) 
are physically implausible on the basis of kinematical parameters extracted 
from galactic orbits and average abundances for the different PN types. 

Recent work on the much better known mass distribution of white dwarfs 
fully support this concl usion, as can be seen from the white dwarf mass distri- 
bution of Madej et al. (|2004 ) , especially regarding the largest mass observed 
in the sample, which in our c ase is under 0.8 Mq. According to the investi- 
gation by Madej et al. (|2004l ), which was based on a large sample of about 
1200 white dwarfs from the Sloan Digital Sky Survey, their mass distribution 
peaks around 0.56 Mq, with very few objects having masses larger than about 
0.8 Mq, in excellent agreement with our results. According to this work, less 
than 5% of the central stars are expected to have masses larger than 0.8 Mq, 
a muc h lower fraction than o bserve d either in the samples by Kudritzki et al. 
|l997h and Pauldrach et al. |20o4 ). 

The mass distributions of CSPN and white dwarfs have also been recently 
considered by Gesicki & Zijlstra (|2007h . based on a dynamical method which 
allows mass determinations within 0.2 Mq. It results that both the CSPN 
and white dwarf distri butio ns peak around 0.6 M as in Madej et al. (|2004j ) 
and Stasinska et al. (| 19971 ) . although the white dwarf distribution shows a 
broader mass range. The CSPN distribution has essentially no objects with 
masses higher than 0.7 M in their sample, while for the DA white dwarfs - 
presumably the progeny of H-rich CSPN - a small fraction of object s hav e 
masses larger than 0.8 Mq, similar to the results of Madej et al. (|2004j ) . 
Although there may be some differences between the recently obtai ned white 
dwarf mass distributions, as discussed by Gesicki & Zijlstra (2007J), they all 
agree in the sense that any sizable sample of CSPN is expected to have masses 
close to 0.60 M Q , a result that is reproduced by our N/O masses. 

A fur ther s upport to the N/O mass es com es from the revision by Tinkler & 
Lamers (|2002l ) of the Kudritzki et al. (jl997l ) masses based on a homogeneous 
set of parameters derived from Zanstra temperatures, dynamical ages of the 
nebulae and evolutionary tracks. According to these results, an average mass 
of 0.60 M Q is obtained, in excellent agreement with the results of Table 1, and 
showing no objects with masses larger than 0.8 Mq. 

Finally, it could be mentioned that 5 central stars of PN in the galactic 
bulge have been recently studied by Hultzsch et al. (|2007l) based on high- 
resolution Keck spectra and NLTE modelling. In view of the fact that the 
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distance to the galactic bulge is well known, the stellar parameters can be 
determined with better accuracy than in the case of field PN, although there is 
some uncertainty due to the adopted extinction law. The masses derived from 
their "Method 1" , which seems to be more reliable, come from evolutionary 
tracks on the logT e // x logg diagram, and are lower than about 0.80 Mq in 
all cases. Also, the estimated luminosities are close to the values obtained by 
Eq. (3), supporting our results given in Table 1. 



3. THE MODIFIED WIND MOMENTUM OF CSPN 

For a star with radius i?* having a wind with terminal velocity Voc and a 
mass loss rate dM/dt, the modified wind momentum is defined as 



Pu 



dM 
~dT 



(4) 



It is well known from radiation driven wind theory that the mass loss rate 
is a function of the stellar luminosity L„ and the effective mass M e f / , which 
can be written as 



dM 1/Q s=k 

~d~t * M « t 



(5) 



where a is the so-called force multiplier parame ter, w hich is about a ~ 2/3 fo r 
hot stars (see for example Lamers & Cassinelli& Pauldrach et al. Q. 
The exponent a is the power law exponent of the line strength distribution 
function. The effective mass M e f / is simply the stellar mass M* corrected for 
the electron scattering radiative force, M e // = M* (1 — r e ), where T e is the 
electron scattering efficiency. 

As a first approximation, can be assumed to be proportional to the 
stellar escape velocity v esc = (2GM e ///i?*) 1 ' 2 , so that Woo oc y/M e ff/R*. 
For exam ple, fo r a large sample of ga lactic CS PN f rom a variety of source s 
fPerinotto [l99l Kudritzki et al. Il997t Malkov fl997t Pauldrach et al. |2004) . 
we obtain a scaling factor between the terminal and escape velocity of about 
2 to 4. Combining this equation with eqs. (4) and (5), it is easy to see that 
the modified wind momentum depends weakly on the stellar mass, so that 
a relation o f the form p w oc L* is obtained. For instance, Pauldrach et al. 



(|2003l . 120041) obtained a relation that produces a very good fit to the properties 
of hot massive galactic stars given by 



log Pw = 20.479 + 1.507 log(L*/L c . 



(6) 



where p w is in g cm s~ 2 . Fig. 3 show s the modified wind momentum with 
CSPN data by Kudritzki et al. (|1997l ) (empty circles) and Pauldrach et al. 
(|2004f ) (filled circles), using the same symbols as before, as a function of the 
stellar luminosity . The solid line represents the fit to hot, massive stars by 
Pauldrach et al. 1)20041 ) given by Eq. 6. These stars are not shown in Fig. 3, 
and occupy the upper right corner of the figure. Note that for NGC 2392, 



METALLICITY EFFECTS IN CSPN WINDS 



9 



en 

CD 
O 

s 

Cl 

CT> 
O 



o 



CO 
CN 



CO 



o Kudritzki A Kudritzki corrected 
• Pau drach a Pou droch corrected 








yr 




A ° ^ ' 






4a a> -v • Z^^^ 






^ - " o . ^ " 







log ( 



Fig. 3. The modified wind momentum as a function of the stellar luminosity. Sym- 
bols are as in the previous figures: e mpty c ircles: data by Kudritzki et al. (1997); 
filled circles: data by Pauldrach et al. (|2004T ); emtpy and filled triangles are the cor- 
responding corrected values adopting the N/O masses. The metallicity dependence 
of the modified wind momentum is illustrated by the dashed lines, corresponding to 
e(O) = 7.5 and 9.2 (see text). 



NGC 3242, IC 4637 and Tc 1 the mass loss rates from Kudritzki et al. (|1997fl 
are upper limits. 

Adopting the new N/O stellar masses, as shown in the previous section, 
and th e luminosities from the mass-luminosity calibration by Kudritzki et al. 
(|1997l ). the stellar radii can be obtained, keeping the effective temperatures 
originally determined. Therefore, the modified wind momentum p w can be 
determined using the originally derived M and Voc and plotted as a function 
of the luminosity. This is also shown in Fig. 3, where the the ope n tria ngles 
correspond to our correction of the results by Kudritzki et al. (119971) and 
filled triangles refer to M and data by Pauldrach et al. (|2004l ). It is 
clear that this is not a strictly consistent procedure, as the main stellar and 
wind parameters (effective temperature, gravity, abundances, radius, terminal 
velocity and mass loss rate) should in principle be derived in a consistent way. 
However, most of these parameters are essentially unchanged in both analyses. 
In fact, a direct comparison of the Kudritzki et al. and Pauldrach et al. results 
shows that there is a general agreement for logg, Uqo, dM/dt, and T e ff. 
Only the stellar mass and luminosity show larger discrepancies: the masses 
by Pauldrach et al. are generally larger than those by Kudritzki et al., while 
the opposite is true for the luminosities. 

From Fig. 1 we can see that the new N/O masses of CSPN have a more 
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restricted range, so that the same can be expected for the luminosities, as 
confirmed by Fig. 3. From Fig. 3 we can see also that the modified wind 
momentum has now a less tigh t corr elation with the luminosity, similar to 
the results by Kudri tzki e t al. (|1997l) and in disagreement with the results 
by Pauldrach et al. (|2003l ; 12004 ). It can also be seen that the same results 
are obtained for both sources of data, that is, the empty and filled triangles 
in Fig. 3 show essentially the same dispersion. As mentio ned i n Sect. 2, 
these results are in agreement with the results by Napiwotzki (2006), who has 
argued against the large masses for CSPN found by Pauldrach and co-workers, 
on the basis of the use of kinematic parameters obtained from galactic orbits 
and average abundances of the different PN types. 

It should be mentioned that the procedure adopted here does not signif- 
icantly affect the stellar gravity logg, which is a parameter generally well 
determined from the observations. In fact, a detailed comparison of the re- 
sulting logg values with t hose origina lly derived either by Kudritzki et al. 
(1997) or Pauldrach et al. ( 20031 : 2004 ) shows an excellent agreement, in the 
sense that (i) our values lie between the original values determined by these 
sources, and (ii) the a verag e deviations in log g am ount to 0.2 a nd 0.1 dex for 
the Kudritzki et al. |l997h and Pauldrach et al. |2003t |2004) data, respec- 
tively. An alternative approach would be to derive the masses from the N/O 
abundances and the stellar radius from the original logg values, so that the 
luminosity would follow from the observed effective temperature. However, we 
preferred the procedure described earlier, since the mass-luminosity relation 
is well founded theoretically, and any modifications in the adopted calibration 
could be easily introduced in this procedure. 

The N/O masses depend on the measured relative abundances of nitro- 
gen and oxygen, so that uncertainties in these measurements would affect the 
derived masses and consequently, the luminosities. However, N and O abun- 
dances relative to hydrogen are among the best determined abundances in 
photoionized nebulae, with uncertainties generally lower than 0.2 dex. Abun- 
dances relative to oxygen, such as the N/O ratio needed by our procedure, 
are generally better determined, usually within 0.1 dex. Adopting this as an 
average, Eqs. (1) and (2) would lead to a maximum uncertainty of about 
10% for central star masses in the range 0.6 < M/Mq < 0.8. Clearly, such 
an uncertainty would have a minor influence in the derived luminosities as 
plotted in Fig. 3. 



4. METALLICITY DEPENDENCE OF THE MODIFIED WIND 
MOMENTUM OF CSPN 

4.1. The modified wind momentum 

A key point in the analysis of the modified wind momentum is the ex- 
pected dispersion in the plo t of lo g p w x logL/L Q for a given luminosity. The 



results by Kudritzki et al. (|1997f ) show a considerable dispersion for CSPN, 
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while the more recent results by Pauldrach et al. |20cd |2004() . and Paul- 



drach (|2005l ) obtain a much lower dispersion and, consequently, a better fit 
to the line associated with the hot, luminous stars. Our results using the 
N/O masses indicate a larger dispersion in log p w for both sources of data 
on the stellar winds, which is partly a consequence of our restricted mass - 
and therefore luminosity - interval. It is then interesting to investigate the 
expected dispersion at a given luminosity, which we will assume to be caused 
by the dispersion in the stellar (and nebular) metallicity. 

The metallicity dependence of win ds fro m red supergiants and AGB stars 



was recently reviewed by van Loon (|2006f ). The mass ejection process in 



these stars is assumed to be t he act ion of the stellar radiation pressure on 



grains. According to van Loon (2006), the dust to gas ratio tp is proportional 
to the stellar metallicity, ip oc Z. Since the grain absorption optical depth 
t is proportional to the dust to gas ratio, r oc ip, we have r oc Z. It is 
reasonable to assume the same is true for CSPN, here attributing the optical 
depth to the absorption in ressonance metallic lines, so that we have that 
p w oc Mvoa oc r oc Z k , where we may assume for simplicity k ~ 1, combining 
any optical depth effects th at eith er dM/dt or might have. In the 
case of red giants, van Loon (|2006T ) attributes the metallicity dependence on 
the expansion velocity, leaving the mass loss rate essentially independent of 
the metallicity. For CSPN, however, this may not apply, since the radiation 
pressure acts directly on the absorption lines, and is not a continuum-opacity 
process as in red supergiants and AGB stars. 

A similar relation, namely, p w oc Z 9 has also bee n sugg ested for galactic 



hot stars with log L/L© < 5.7 by Lamers & Cassinelli (|1996f) as an upper limit 
for the correction of the mass loss rates due to metallicity. Such a relation is 
similar to the metallic ity de pendence of Wolf-Rayet winds, as recently studied 



by Vink & de Koter (|2005h on the basis of Monte Carlo calculations. As an 



example, for the less evolved WN stars in which the wind is due to the action 
of the radiation pressure on many metallic lines, the metallicity dependence 
is similar as for galactic O stars, namely, dM/dt oc Z m , where m ~ 0.86. For 
WC winds they find m ~ 0.66. 

Adopting p w oc (i„/L Q ) a , as discussed in Section 3, where a — 1/a, we 
can introduce a dependence on the metallicity Z by writing 

p w oc Z ( y— ) • (7) 



L 



In the case of CSPN, the wind strength is basically related to the abundances 
of highly ionized heavy elements, which produce P Cyg profiles of optical and 
ultraviolet lines. Iron is the usual metallicity indicator in stars, but the iron 
lines are generally weak and the Fe abundance is depleted in photoionized and 
planet ary ne bulae due to grain_ formation (cf. Perinotto et al. 1999; Dcctjcn 
et al. 1999t Pottasch et al. 20081 ). so that the abundances of this element 



measured in the nebulae are usually taken as lower limits. However, the 
metallicity of PN can be accurately determined on the basis of measurements 
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of the oxygen abundance, as well as of other elements such as Ne, S, etc 
all of which are strongly correlated, as shown by many recent investigations 
both in the Milky Way and in other galaxies (cf. Richer & McCall 



200 



Maciel et al. 120061 ). As recently discussed by Idiart et al. (|20071 ), average 



6; 



[O/Fe] x [Fe/H] relationships can be obtained from stellar data or theoretical 
models, which may provide direct O(Fe) relations. Also, recent work on the 
[Fe/H] radial gradient in the galac tic di sk obtained from open cluster stars 
and cepheid variables (Maciel et al. 2005) shows a strong correlation with the 
gradients of O, S, etc. determined from photoionized nebulae, thus confirming 
these elements as representative of the nebular metallicity. Therefore, we can 
safely assume that in PN and hydrogen-rich CSPN the metal abundance Z is 
proportional to the oxygen abundance by number O/H, that is, Z = k (O/H). 
For example, for the solar photosphe re we have e(O) = log(O/H) + 12 ~ 



8.66 and Z ~ 0.012 (Grevesse et al. 120071) . so that k ~ 26. Therefore, the 



modified wind momentum can be written as 



where (3 is a constant to be determined. We have then 

log p w = log + a log (j^j + lo S 
Using the standard definition e(O) = log(0/H) + 12, we have 



(8) 



(9) 



\ogp w = log (3- 



[e(O) - 12] 



(10) 



According to Mendez (|1991r ). all CSPN considered in this work are H-rich 
objects, in which stellar H features are clearly distinguished in the CSPN 
spectrum. Their H and He abundances are normal, with the possible excep- 
tion of NGC 2392, which shows some He excess. All the other objects have 
"normal" helium a bundances of about 10% by number (for the actual values 
see Mendez et al. 1988 . 1990l) . again suggesting that these stars have normal 
Fe abundances, and that their metal abundances are well correlated with the 
oxygen abundance. Most of these stars are of spectral type Of(H), so that the 
He II 4686 line appears as a narrow emission line, and two are of type O(H), 
in which He II 4686 appear in absorption. In the sample of 115 CSPN studied 
by Mendez |l99ll ). 62% are H-rich and 33% are H-deficient. These objects 
have probably left the AGB with a H-rich photosphere, and are on th eir wa y 
to become DA white dwarfs. As can be seen from figure 2 of Mendez (|199ll ). 
most H-rich CSPN have "normal" He abundances, that is, He/H ~ 0.10, so 
that their chemical composition includes about 90% of H by number, as in 
the case of the nebulae themselves. Therefore, we can safely adopt average 
oxygen abundances in these objects as in the nebulae themselves, that is 



e(O) = e(PN) = e(CSPN) ~ 8.50 , 



(11) 
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TABLE 2 

EXPECTED DISPERSION OF THE MODIFIED WIND MOMENTUM. 



e 2 (0) 


A log p w 


e a (0) 


A log p w 


9.20 


0.70 


8.30 


-0.20 


9.10 


0.60 


8.20 


-0.30 


9.00 


0.50 


8.10 


-0.40 


8.90 


0.40 


8.00 


-0.50 


8.80 


0.30 


7.90 


-0.60 


8.70 


0.20 


7.80 


-0.70 


8.60 


0.10 


7.70 


-0.80 


8.50 


0.00 


7.60 


-0.90 


8.40 


-0.10 


7.50 


-1.00 



which is essentially the same average oxyg en abundances of galact ic HII re 



gions (see for example Costa & Maciel l2006[ Henry & Wort hev 1 19991 ). We can 



identify Eqs. 10 and 6, from which we have a — 1.507, and log (3+[e(0) — 12] = 
20.479, so that log (3 = 23.98. Therefore we have 

logp„ = 23.98 + 1.507 log (j^j + [e(O) - 12] , (12) 

again with p w in cgs units (g cm s -2 ). The new p w values are shown in the last 
column of Table 1. For a given luminosity, the expected dispersion in log p w 
due to a difference in the metallicity is simply Alogp^ = Ae(O). The two 
dashed lines in Fig. 3 show the results for e(O) = 7.50 and e(O) = 9.20, which 
correspond approximately to the minimum and maximum observed values of 
the oxygen abundance, respectively. From the metallicity distribution of PN 
and HII regions (as well as cepheids and CSPN), we know that the large 
majority of objects have 8.00 < e(O) < 9.0, so that the derived dispersion 
is very well defined. Table 2 shows the dispersion A log Piu &-S 0j function of 
the adopted value of the oxygen abundance e(O). It can be seen that, for the 
average e(O) = 8.50 and 7.5 < e(O) < 9.2, the spread is from Alogp,,, = 
— 1.00 to 0.70, that is, a total spread of 1.70 can be expected, which is similar 
to the observed scatter in Fig. 3 using the N/O masses and both the Kudritzki 
or Pauldrach databases. 



4.2. Empirical evidences 

Empirical evidences of metallicity effects on the modified wind momentum 
of CSPN are difficult to obtain for several reasons. First, the modified mo- 
mentum probably shows a stronger dependence on the stellar luminosity, so 
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that any metallicity variation is expected to become apparent only for stars 
having similar luminosities. Second, the actual stellar luminosity has some in- 
trinsic uncertainty, which basically derives from uncertainties in the distance 
scale, which affects particularly planetary nebulae and their central stars. Fi- 
nally, the abundances themselves have an uncertainty of typically 0.2 dex, 
depending on the element considered. Despite these shortcomings, data on 
some of the objects analyzed in this paper show some hints in the sense that 
an enhanced metallicity is correlated with a higher wind momentum. 

Taking into account the results shown in Table 1, we notice that the objects 
He2-108 and He2-131 have approximately the same luminosity, logL/L Q ~ 
3.95. The first object is clearly more metal-poor than the second, with a 
difference of abou t 0.5 dex, or almost a factor 3. According to the results of 



Kudritzki et al. (|19970 . which we favour here, He2-131 has a much higher 
modified wind momentum (logp w ~ 27.82, cgs units) than He2-108 (logp w ~ 
27.28), which may be interpreted as a metallicity effect. Using our "corrected" 
values given in Table 1 the same pattern is observed. Also for the pair IC 4593 
and Tc 1, Table 1 shows similar luminosities, namely, logL/L© ~ 3.77. Again, 
the first object is more metal-poor than the second by a factor of about 2. This 
pattern is repeat ed in the determination of the modified wind momentum by 



Kudritzki et al. (j 19971 ) and in our c orrected values given in Table 1, although 



the value by Kudritzki et al. (jl997l ) for Tc 1 is strictly an upper limit. 



If we consider now the objects studied by Perinotto (|1993f ). we notice that 
NGC 1535 and NGC 7009 have similar luminosities, logL/L Q ~ 4.0, but the 
former is more metal-poor than the latter by a factor of 1.7. Again the same 
tenden cy is repeated in the modified wind momentum, for which Perinotto 



(1993) obtains logp lu ~ 25.7 and 26.0, respectively. 



An indirect empir i cal ey idence can also be obtained from the sample stud- 



ied by Malkov (jl997l : [1998) , in a self consistent determination of several pa- 
rameters of galactic planetary ne bulae and their central stars. Among other 
parameters, Malkov dl997i : Il998h obtains central star luminosities and neb- 



ular abundances for about 130 nebulae in the galactic disk. If we make the 
reasonable assumption that the modified wind momentum p w of these objects 
depends on the stellar luminosity in the same way as the small sample of 
CSPN with well studied winds, we could use an equation such as eq. (6) to 
estimate this quantity. We notice that 74% of the galactic disk nebulae in the 
sample have oxygen abundances higher than the average value adopted here, 
e(O) = 8.50, while 26% of the nebulae are metal-poor, having e(O) < 8.50. 
Estimating the modified wind momentum from eq. (6), we obtain similar fig- 
ures, that is, 77% of the stars have \ogp w > 25.0 in cgs units, while 23% 
have lower values of this quantity, namely, logp^ < 25.0. Of course, this 
is an indirect evidence that may be affected by the assumption made and 
the uncertaintes involved in the determination of these quantities, but it is 
reassuring that some correlation between the modified momentum and the 
nebular metallicity is suggested. 
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5. SUMMARY AND CONCLUSIONS 

In this work we have considered a sample of well-studied CSPN for which 
the wind properties as well as the basic stellar properties are known, and 
analyzed two different problems. 

First, we have confirmed that there is a large discrepancy between the 
stellar masses as derived from standard mass-luminosity relation for AGB 
stars and those derived from stellar atmosphere models affected by winds. 
We have then determined new values for the central star masses based on a 
correlation of the core stellar mass and the nebular N/O abundances, and 
determined the masses of the central stars. We showed that the masses agree 
very well with the independent estimates of the masses of CSPN, and that 
the observed mass spread is small, entirely consistent with the expected mass 
distributions of CSPN and white dwarfs. 

Second, we have used the newly derived N/O masses and the standard 
mass-luminosity relation in order to investigate the dispersion of the modi- 
fied wind momentum at a given luminosity. We have taken into account two 
different sets of stella r properties derived from different model atmospheres 



(Kudritzki et al. 119971 and Pauldrach et al. 120041 ) , and showed that the disper- 
sion on the p w x plane is real, and can be well explained by the metallicity 
dispersion observed in the nebulae, which are a reflection of the properties 
of the central stars. Similar results can al so be obtained by considering the 
set of objects analyzed by Perinotto |l993h . in which case the mass loss rates 
have been determined with the SEI (Sobolev plus exact integration) and CFS 
(comoving frame) method. 
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